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ABSTRACT 

Matter accreting onto black holes has long been known to have standing 
or oscillating shock waves. The post-shock matter puffs up in the form of a 
torus, which intercepts soft photons from the outer Keplerian disc and inverse 
Comptonizes to produce hard photons. The post-shock region also produces 
jets. We study the interaction of both hard photons and soft photons, with on- 
axis electron-positron jets. We show that the radiation from post-shock torus 
accelerates the flow to relativistic velocities, while that from the Keplerian 
disc has marginal effect. We also show that, the velocity at infinity or terminal 
velocity depends on the shock location in the disc. 

Key vifords: Accretion, accretion discs - black hole physics - radiation mech- 
anism: general - radiative transfer - ISM: jets and outflows 

1 INTRODUCTION 

Jets around quasars and micro-quasars show relativistic terminal speeds. While jets are quite 
ubiquitous and are associated with a wide range of celestial objects, only some jets around 
quasars and micro-quasars show highly relativistic terminal speed {e.g., GRS 1915-fl05, 
Mirabel & Rodriguez 1994; 3C 273, 3C 345, Zensus et al. 1995; M87, Biretta 1993). These 
relativistic jets are generally associated with compact objects and circumstantial evidences 
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show that many of these central gravitating objects are black holes. Black holes do not have 
'hard surfaces' nor do they have atmospheres, hence if observations show that many of these 
jets come from the vicinity of the black hole, then they must originate from the accretion 
discs around these black holes. 

Inner boundary conditions of matter accreting onto a black hole are (i) supersonic and 
(ii) sub-Keplerian. Liang and Thompson (1980) showed that sub-Keplerian matter accreting 
onto a black hole has at least two X-type critical points. In much of the parameter space, 
it has been shown that the supersonic matter crossing the outer critical point, under goes 
centrifugal pressure mediated shock (Fukue, 1987; Chakrabarti, 1989), becomes subsonic 
and enters the black hole through the inner X-type critical point. Entropy is generated at 
the shock making the post shock region hot. The region in which the flow slows down may 
be extended if the shock conditions are not satisfied. This hot, slowed down region is puffed 
up in the form of a torus (hereafter, CENBOL = CENtrifugal pressure supported BOund- 
ary Layer; see Chakrabarti et al. 1996, hereafter CTKE96). This disc due to its advection 
term may be called advective accretion disc (Chakrabarti 1989, hereafter C89; Chakrabarti 
1990, hereafter C90; Chakrabarti 1996, hereafter C96; CTKE96). Chakrabarti and Titarchuk 
(1995, hereafter CT95), proposed a disc model which contains both the Keplerian matter 
and the sub-Keplerian matter. In this model, the Keplerian matter is of higher angular 
momentum and low specific energy, and settles around the equatorial plane to form the Ke- 
plerian disc (see, Shakura & Sunyaev 1973; Novikov & Thorne 1973, hereafter NT73) while 
the sub-Keplerian matter with high energy and lower viscosity flanks the cooler Keplerian 
disc from the top and bottom, sandwiching the Keplerian disc known as the sub-Keplerian 
halo (see, CT95; Chakrabarti 1997, hereafter C97) in the literature. The sub-Keplerian halo 
may suffer a standing shock at x^, a few tens of Schwarzschild radii, and it may be sustained 
there if the post-shock thermal pressure is sufficiently high. The shock compresses the pre- 
shock flow making it denser and at the same time hotter. In the model solution proposed 
by Chakrabarti & Titarchuk (CT95), the post-shock region is comprised of a mixture of the 
Keplerian and sub-Keplerian components. Thus, though the sub-Keplerian halo (pre-shock) 
is optically thin for the radiations from Keplerian disc, the post-shock torus could be opti- 
cally thin, intermediate, or even thick depending on the Keplerian and sub-Keplerian rates 
(see CT95, CTKE 96, C97 for details). This is because (1) the Keplerian radiation falls 
on it at a glancing angle, thereby increasing the path length and (2) the mixed matter in 
this region has higher density. CT95 showed soft radiation from the cool Keplerian Disc are 
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Figure 1. Cross-sectional view of Two Component Accretion Disc Model. Only the top half is shown. 

inverse- Comptonized by the CENBOL to produce the hard radiation. If the sub-Keplerian 
halo rate {Mh) is higher, then it supphes more hot electrons to the CENBOL than the soft 
photons from the Keplerian disc and hence the soft photons cannot cool down the CENBOL 
significantly. Thus CENBOL remains puffed up and hot and can intercept a large number of 
soft photons and inverse- Comptonize them to produce the hard power-law tail of the accre- 
tion disc spectrum — a state called hard state. If, on the other hand, the Keplerian accretion 
rate (M^) is higher, then it supplies more soft photons to cool down the CENBOL region. 
This results in more power to the soft end of the accretion disc spectrum — a state known 
as the soft state. Recently, however, Chakrabarti & Mandal (2003) showed that raising the 
Keplerian rate even higher not necessarily softens the spectrum, for, the Keplerian flow also 
adds to the number density of electrons in the post-shock region and at some point, the 
spectrum starts to be hardened once more. 

This kind of hybrid disc structure is known as the Two Component Accretion Flow or 
the TCAF disc (see, CT95, CTKE 96, Ebisawa et al. 1996, C97), and has wide observational 
support (Smith et al. 2001, 2002). In Fig. (1), a schematic diagram of such a disc is presented. 
The Figure shows the Keplerian disc, is sandwiched by the sub-Keplerian halo. The shock 
location {xg) and the compact object are also shown. Jets are shown close to the axis of 
symmetry. Thus jets are illuminated by the Keplerian disc with soft photons and by the 
CENBOL with hard photons. 

Chakrabarti and his collaborators have also showed that the CENBOL can drive a part of 
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the infalling matter along the axis of symmetry to form jets (Chakrabarti, 1998; Chakrabarti, 
1999; Das & Chakrabarti, 1999; Das et al. 2001). There are wide support that the jets are 
indeed coming out from a region within 50 — 100 Schwarzschild radius of the black hole 
(Junor, Biretta & Livio, 1999). Similarly, it is beheved that jets are produced only in hard 
states (see, Gallo, Fender & Pooley, 2003, and references therein). Thus it is natural to study 
interaction of hard radiation from the CENBOL and the outflowing jets, with the particular 
interest of studying, whether momentum deposited to the jet material by these hard photons 
can accelerate them to ultra-relativistic speeds. 

Investigation of interaction of radiation and astrophysical flow is not new. A number of 
astrophysicists have directed their efforts in this particular fleld of study while the consider- 
ation of the associated accretion disc depended on their personal choice or the popularity of 
the particular model of accretion at the given time. Icke (1980) studied the effect of radiative 
acceleration of the gas flow above a Keplerian disc. But the effect of radiation pressure on 
the gas flow was ignored. Sikora and Wilson (1981) showed that even if the jet is collimated 
by geometrically thick discs (Lynden-Bell, 1978; Abramowicz and Piran, 1980), radiation 
drag is important for astrophysical jets. Piran (1982) while calculating the radiative accel- 
eration of outflows about the rotation axis of thick accretion discs, found out that in order 
to accelerate outflows to 7 > 1.5 (where 7 is the bulk Lorentz factor), the funnels must be 
short and steep, but such funnels are found to be unstable. Sol et. al. (1989) proposed a 
two-flow model for jets, one consists of relativistic particles (electrons and positrons) and 
of relativistic Lorentz factor, while the other being normal, mildly relativistic plasma. In 
a very important paper, Icke (1989) considered blobby jets about the axis of symmetry of 
thin discs and he obtained the 'magic speed' of — 0.451c (c — the velocity of hght), 
Vmheing the upper limit of terminal speed. Fukue (1996) extended this study for rotating 
flow above a thin disc and drew similar conclusions, although for rotating flow, away from 
axis of symmetry the terminal speed was found out to be a little less than the magic speed 
of Icke. To summarize, earlier works showed that it is not possible to accelerate jets to ultra 
relativistic terminal velocities, by radiations from the earlier accretion disc models. What is 
more discouraging is the existence of moderate levels of equilibrium speed {veg i.e., speeds 
above which there would be radiative deceleration). Recently Fukue et al. (2001) has done 
similar investigations on interaction of radiation and pair plasma jets, but they took a disc 
model which consisted of inner ADAF region (non luminous) and outer slim disc (luminous), 
which resulted in a relativistic terminal speed. 
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We are working in a different regime, i.e., less luminous Keplerian disc and more luminous 
post-shock torus or CENBOL. Since hard radiations are expected to emerge out of the 
optically thin CENBOL, its intensity (photon counts per unit area per unit time) is low, but 
it 'looks' directly into the jet vertically above and hence eventually deposit its momentum 
into the later. Radiation from a hot CENBOL is likely a source of pair-production and hence 
the possibility of radiative momentum deposition is likely to be higher even for radiations 
from CENBOL hitting the outflow at an angle (sec, e.g., Yamasaki, Takahara, and Kusunose, 
1999 for the mechanism of pair-production from hot accretion flows.). This is why we believe 
that the direct deposition of momentum may be important. Chattopadhyay & Chakrabarti 
(2002a) showed that hard radiations from the post-shock region (CENBOL) do accelerate 
electron-proton plasma to mildly relativistic terminal speeds. Chattopadhyay & Chakrabarti 
(2002b) also reported that hard radiations from the CENBOL do not impose any upper limit 
for terminal speed. 

In the present paper, we solve the equations of photo-hydrodynamics of jets for radiations 
coming out from the TCAF discs, where the radiation fields from both the inner CENBOL 
and the outer Keplerian thin disc are considered. We show that, while the equilibrium 
velocity closer to the black hole depends on Xg, and the ratio between Keplerian disc and 
CENBOL luminosities, the terminal speed or the jet velocity at infinity depends on the 
relative proportions and also on the actual magnitude of various moments of radiation. 
We also show that, in hard states (in our parlance, CENBOL radiation dominating over 
Keplerian radiation), optically thin jets can be accelerated to ultra-relativistic speed. 

In the next Section, we present the model assumptions and the equations of motion and 
compute various moments of radiation field. In §3, we present our solutions and finally in 
§4, we draw our conclusions. 

2 ASSUMPTIONS, GOVERNING EQUATIONS AND COMPUTATION OF 
THE MOMENTS OF RADIATION FIELD 

2.1 Assumptions and Governing Equations 

In our analysis, the curvature effects due to the presence of the central black hole mass is 
ignored. The metric is given by, ds'^ — c^dt^—dr'^—r'^d(t)^—dz^, where, r, 0, and z are the usual 
coordinates in cylindrical geometry and ds is the metric in four-space. The four- velocities are 
u^. We follow the convention where the Greek indices signify all four components and the 
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Latin indices represent only the spatial ones. The black hole is assumed to be non-rotating 
and hence the strong gravity is taken care of by the so-called Paczyhski-Wiita potential 
{e.g., Paczynski & Wiita, 1980). 

We also do not consider generation mechanism of jets. As the astrophysical jets are 
observed to be extremely coUimated (Bridle & Perley, 1984), and generally aligned along 
the normal to the host galaxy, we assume the jet to be along the axis of symmetry. Thus 
the transverse structure of the jet is ignored, i.e., u'^ = u'^ = and d/dr = d/d(f) = 0, where 

and are radial and azimuthal components of four velocity. We are looking for steady 
state solutions. Hence, d/dt — 0. We also assume the gas pressure is neghgible compared 
to the radiation pressure. This is perhaps the case especially inside the funnel wall close to 
the axis. The derivation of the equations of motion of radiation hydrodynamics for optically 
thin plasma, was investigated by a number of workers. A detailed account of such derivation 
has been presented by Mihalas & Mihalas (1984; hereafter MM84) and Kato et al. (1998; 
hereafter K98), and are not presented here. Enforcing the above assumptions, the equation 
of motion presented in MM84 and K98, reduces to; 

= -7^^ + — [7— - -f^Eu' - u'P'' + uH2^u'— - u'u'P'')] , 

az [z — Zy m V c c J 

where, is the 2;-component of four velocity, G, Mb, ctt, and m are the universal gravitation 
constant, mass of the black hole, Thomson scattering cross-section and mass of the gas 
particle, respectively. E, and P^^ are the radiative energy density, radiative flux and 
radiative pressure on the axis of symmetry, and 7(7 = Ut) is the Lorentz factor. The above 
equation can be re-written as, 

^2 + - - + «'(27«'^ - u'u'V)] , (1) 



du' GMb 



dz {z - 2f 
where, £ = ^E, T = and V = ^P^\ 

For simphcity, we will not compute the shock location Xg or the the CENBOL luminosity 

{Lc) - instead, we will supply them as free parameters. They can be easily computed from 

accretion parameters {e.g., C89, CT95, Das et al. 2001, Chattopadhyay et al. 2003). We 
assume that the outflow is made up of purely electron-positron pair plasma. 

2.2 Computation of radiative moments from TCAF disc 

The radiation reaching each point on the jet axis, is coming from two parts of the disc, 
namely, the CENBOL and Keplerian disc, hence all the radiative moments should have 
both the contributions. 
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Figure 2. Schematic diagram of a two component accretion flow (TCAF). O is the position of the black hole. Centrifugal 
pressure dominated boundary layer (CENBOL) is the puffed up region between O and Xs, the shock location. Thick line 
between Xs and Xo is the Keplerian disc where Xo is the outer boundary of the Keplerian disc. DC and D'C' are the local 
normals of the CENBOL and Keplerian disc. B{0, z) is the field point where various moments of radiation fields are computed. 
D{r', z') is source point on the CENBOL and D'{r'j^,0) is the source point on the Keplerian disc, r' is radial coordinate of the 
CENBOL and is that of the Keplerian disc. The sub-Keplerian halo is not shown. 



In Fig. (2), a schematic diagram of the cross-section of the disc structure is presented. 
The black hole is situated at O. The region bounded between O and Xg is the CENBOL, 
and the thick line between Xs and Xg is the Keplerian disc, where Xo is the outer boundary 
of the Keplerian disc. Thus Xg is the outer boundary of CENBOL and inner boundary of 
the Keplerian disc. The inner boundary of Keplerian disc, as seen from the position B is 
Xxi- The shock height is hs'^asxl/'^{xs — 1), where Ug is the equatorial sound speed at shock 
and depends on Xs- In other words hs — hs{xs), but as we are not solving the accretion disc 
equations simultaneously, we have to make some estimate of hs which will closely mimic 
reahty. Chakrabarti solutions (C89, C90, CT95, C96) show that if Xg — lOr^ (where rg is 
the Schwarzschild radius), then the temperature at the shock is Tio~1.56xlO-'^-^X. Assuming 
that the shock temperature to be Tg = Tio(10/a;s), one can estimate the shock height to be 
hs'^0.6{xs — 1). The inner surface of the CENBOL is assumed to be conical. The radiations 
from a point D{r', z') on the surface of CENBOL is primarily along the local normal DC. 
Similarly, the radiation from a point D'{r'j^, 0) on the Keplerian disc, is along its local normal 



D'C 



The radiative moments at B are; 
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— = - f ivdn = - ( / icVcd^c + / iKiUnK) , (26) 

and 

P^o ^ 1 y ifVd^ = ^ [J JcV-cVc^^c + J ^iKiyidflK^ . (2c) 

In the Eqs. (2a-2c), /, dJl, are the frequency integrated intensity from the disc, differ- 
ential sohd angle at B and T are the direction cosines at B, for example = {z — z')/BD 
and = z/BD' . Suffix C and K represent quantities linked to CENBOL and the Keplerian 
disc respectively. The expressions of solid angles subtended at B from D and D' are given 
respectively by, 



and 



d^c = ^2—7 ^cos{lCDB), 3a 

r + [z — z'Y 

dn^ ^ !J<P<^cos{IC'D'B), (36) 

where, 9 is the semi-vertical angle of the inner surface of the CENBOL which, for simplicity, 
is assumed to be constant. It is to be noted that, from Fig. (2) and Eqs. (3a-3b), in general, 

> f*^^ ^^i^ differential area on the disc, dQc > dQx- This implies that 

the contribution from the CENBOL to the total radiation field moment is much more, 
compared to that from the Keplerian disc. This will be clear when the comparative study 
of the moments are presented later [see Fig. (3)]. 

It is clear from Eq. (2b) that, due to the symmetry about the axis, only non-zero 
pressure tensor components are P*^ and F'' = F"^ = on the axis. Only F^^O. As we 
consider jets on or about the axis of symmetry and the transverse structure is ignored, P^^ 
and P'^'^ do not enter the equation of motion as these components are only coupled with 
u^{— 0) and u'^{— 0), i.e., we have to only compute P^^, which is exactly what is seen in 
Eq. (1). 

From Fig. (2), it is clear that for the outfiowing matter within the funnel, i.e., when z < 
[hgXo/ — Xg)], radiations from the Keplerian disc do not reach the electrons because these 
radiations are intercepted by CENBOL. CENBOL will reprocess these intercepted photons, 
and re-emit them in all directions, especially towards the jets because of especial geometry 
and directiveness of the CENBOL surface. Due to shadowing effect mentioned above, the 
radiation from the Keplerian disc reaching B is from xxi to Xo, where xxi — Xsz/{z — hg), 
with the additional constraints, Xs<XKi<Xo and Xxi > 0. 

As far as the CENBOL properties are concerned, we follow CT95, where an effective 
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Radiatively driven electron-positron jets from two component accretion flows. 9 
temperature was computed for the CENBOL and the radiation intensity was chosen to be 
uniform. That is: Ic = Lc/nA = iLEdd/'^A= constant, where Lc and A are the CENBOL 
luminosity and the surface area of the CENBOL respectively. LEdd is the Eddington lumi- 
nosity and £ is the CENBOL luminosity in units of LEdd- The Keplerian disc intensity per 
unit sohd angle is Ik = ^^sSf"" ~ V^) (^^73). Let us now multiply ar/m with Eqs. 
(2a-2c), and then integrate over the whole disc to obtain the following integrated quantities 
for the moments, 

27rcote{r' + (z - r'cote)tane}r'drf f^o 2TTz[rK^ - V^rK^^^]dr'K 



f^" 27:cotU{r + [z - r cotU)tant)^r drf /"^ 
\(z-r'cot9V + r'^]y^ ^ Jx, 



' [{Z - r'cotey + r'2]3/2 ' J^^^ (^2 + ^^)3/2 

= 8coEc{z, rin, Xs) + £koEk{z, Xg, Xo) 

= £c + £k, (4a) 

where, r^n is the inner boundary of the disc. 

f-^s 2ncot9{r' + {z - r'cot9)tane}r'{z - r'cotd)dri 



Jx 



\{z - r'cot&f + r'2]2 

■■'^o 2TTz[rK^ - ^/3rK^^^]zdr'K 



TcoFc{z, rin, Xs) + TkqFk{z, Xs, Xo) 

:Fc + Tk (46) 



CO 



2'KcotQ{r' + (2; - r'cot6)tan6}r'{z - r'cot6)'^dr/ 

[{z - r'cote)^ + r'2]V2 

27rz[r'K^ - V3rK''^^]z^dr'K 



I^K. {z^ + )5/2 

= VcQpciz, Tin, Xs) + VkqPk{z, Xs, Xo) 

= Vc + Vk (4c) 

Constancy of Iq allows us to have analytical expressions for Eq, Fq and Pq- These were 
computed by Chattopadhyay & Chakrabarti (2000, 2002ab) and are not repeated here. We 
choose the unit of length to be 2GMb/c^ — the Schwarzschild radius (r^,), unit of time to be 
2GMbI(? and Me is the unit of mass. Thus the unit of velocity is c. In such units, constants 
in Eqs. (4a-4c) are, 

L3xlO^«£aT . 



and 
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We have written the Keplcrian accretion rate in units of the Eddington accretion rate M^ddj 
i.e., TTiK — Mx/MEdd- It is to be noted that, two of the disc parameters i.e., r^n — l-Sr^ 
and Xo — lOOOrg are kept constant through out the paper. In Fig. (3), the space dependent 
part of the moments E, F and P are compared for two shock locations. In Figs. 3(a-b), the 
contributions from the CENBOL and the Keplerian disk are plotted for Xg = lOr^ and in 
Figs. 3(c-d) the contributions from the CENBOL and the Keplerian disc are plotted, where 
the shock location is Xg = '^•QTg. We see that clS Xg IS increased, the CENBOL contributions 
increases while the Keplerian contributions decrease. We also see another fact that generally, 
in the hard state, the CENBOL contribution to the total radiative moments dominates over 
Keplerian counterparts. We see from Fig. (3a) and (3c), that within the funnel {i.e., z<hs), 
Ec > Fc > Pc- We also see that, because of the relatively small size of the CENBOL, 
Ec^Fc^Pc for smaller vahics of z [e.g., zr^50.5rg for Fig. (3a) and z^92.5rg for Fig. (3c)], 
while we sec Ek^Fk^Pk at much larger distances [e.g., 2;~900rg for Fig (3b) and 2;~1500rg 
for Fig. (3d)]. The small size of CENBOL ensures that the direction cosines Iq^I for smaller 
z. In contrast, the larger size of the Keplerian disc ensures only for large z. From Figs. 
3(a-d), it appears that as Eq^Ek, Fq^Fk, and Pc^Pr- So, if CENBOL is more luminous, 
then CENBOL contributions would dominate the Keplerian ones even at large z. While £ko 
(= J-'ko = Vko) depends only on mx, £co (= ^co = "Pco)) depends on both (. and the 
CENBOL surface area. From Figs. (3a) and (3c), we see that with the increase of Xs, even 
if the radiative moments from CENBOL, say for example, Eq increases moderately, but 
calculations show that Ecq decreases appreciably with Xg {e.g., Eqq for Xg — lOr^ is about 
four times larger than Eqq for Xg — 20rp). So Keplerian contributions might be comparable to 
that due to the CENBOL at large z, depending on specific combinations of all the parameters 
£, tfiK and Xg, especially for large values of Xg and low We also observe that the shadow 
effect of CENBOL ensures that at z<{hgXo) / {xo — a;^), Ek = Fk = Pr = 0. 

3 RADIATIVE ACCELERATION 

In the rest of the paper, we will use geometrical units defined in the last section, but for 
simplicity we will keep the same symbols representing variables as in Eq. (1). We define a 
three-velocity v such that — —UiU^/utU^ — —UzU^/utU^. Thus — and Uz — —jv. 
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Figure 3. Variations of spatial part of radiation moments, E (solid) F (dashed) and P (long-dashed) with log(z). (a) Con- 
tributions from the CENBOL denoted by sufRx C, and (b) Contributions from the Keplerian disc denoted by suffix K, shock 

location is Xs = lOrg. (c) Contributions from the CENBOL for Xg = 20rg, and (d) Contributions from the Keplerian disc for 
Xs = 20rg. The suffixes have the same meaning as in the previous two Figures. 

Under such considerations Eq. (1) takes the form, 

dv ^ - 2(iV + - - IV'P + 7^(2^'-^ - 

dz~ 7% ■ ^ 

The first term in the right-hand side of Eq. (6) is the gravitational term with dimensionless 

Paczynski-Wiita potential and the term in the square bracket is the radiative acceleration 

term. We notice that the radiative acceleration term depends both on v and the radiative 

moments. The first and fourth term in the square bracket are accelerating terms while the 
second, third and fifth terms are the decelerating terms, collectively known as the radiation 
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drag terms. We also see that there is a 7^ term in the denominator of the r. h. s of Eq. (6). 
This term ensures that, apart from radiation drag, outflowing matter will be slowed down 
as i>— >1. 

3.1 Equilibrium velocity 

Let us now discuss the concept of equilibrium velocity. As a concept, equilibrium velocity 
is not a new one and has been extensively discussed by a number of astrophysicists (see, 
Fukue 2003, and references there in for details). We want to study this issue in the context 
of jets in the radiation field of TCAF discs. It is defined in a manner that at v — v^q the 
square bracket term in Eq. (6) is zero. Thus, for v > Veq there is deceleration and for v < Veq 
there is radiative acceleration. 

Putting the square bracket term in Eq. (6) to zero we have, 

J^V% -{S + V)Veq + ^ = 0. (7) 

This is a quadratic equation of Veq whose solution is. 



ie - 1), (8) 

where C = + V)/2J^. Prom Eqs. (4a-5b), we have, 

^ {Ec + Pc) + ^({Ek + Pk) 

where, ( — 6.6xl0~^^c^^. It is to be noted that ^ depends on Xg as well as on the ratio 
ttik/^-i but not separately liriK and £. In case {TfiKCj/^^^-i '^eq is completely determined 
by Eci Fc and Pq-, and ^ has no effect in determining v^q. We also see that ^ does not 
depend on the mass of the gas particles m. Thus ^ (and thus fgq) is the same for both 
electron-proton plasma as well electron-positron plasma, provided Xg and rriK/^ is the same 
in both the cases. It is clear from Eqs. (8-9), that V, Vpa^l. Therefore if the 

CENBOL contribution dominates then, v^q'^l within few tens of Schwarzschild radii above 
the disc plane. If Keplerian radiation dominates, then the condition is achieved at much 
larger distance, as we shall see later. We also see that no outflow is possible i.e., Veq<0, i.e., 
if 

J^<0. (10) 
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Radiatively driven electron-positron jets from two component accretion flows. 13 
From Fig. (3) we see that very close to the horizon, due to the torus geometry of the 
CENBOL JF < 0, hence very close to the black hole, not only enormous gravitational pull 
but also the radiative force pushes matter inward. 
Icke's magic speed: 

Let us now recover an important result from Eq. (8). If we consider a thin-Keplerian disc 
of infinite size, then P|? = P^^ = P'^' = \Ek and F^^ = ^Ek, and there is no CENBOL in 
this particular case. Under such conditions S = 2T = 3V, putting these in Eq. (8) we have, 

Veg — — ^{4: — a/T) = 0.451= magic speed of Icke ! 

Thus we see that, if the radiation field, i.e., the radiative properties of the disc, can be 

prescribed, then v^q can be found out easily. 
Equilibrium speed from a TCAF disc: 

Let us now concentrate on the radiations from TCAF discs. As discussed in the introduction, 
jets are produced in the hard state of the accretion disc when the CENBOL is hotter and 
hard state means more power is on the high energy end of the spectrum. That is, Lc > Lk- 
The Keplerian disc luminosity is given by Lk — r'^fl°2TTlK'27rrKdrK, which is a function of 
TTiK, Xg and Xo, which when integrated can be expressed as. 



2nlK'2'KrKdrK = -rrfiK 

Xs 4 



1 2 3 
+ 



Lem. (11) 



Tk ^Tk V 

Thus the Keplerian disk luminosity in units of the Eddington luminosity can be defined as 
^K = Lx/LEdd- As has been stated before, presently we do not compute Lq. but supply 
it. Typical vahics of £ that we shall employ should, in general, depend on l^: itself because 
i. ~ A£x, where A is the enhancement factor by which incident photon intensity is increased 
due to Comptonization and has a value of around 20 — 30 in hard states (CT95). Thus, for 
instance, if Ik ~ 0.05, a typical value of £ ~ 0.1 — 0.15. However, while we shall choose I 
and £k in these regions, we shall use them as free parameters, since presently we are not 
interested in computing the terminal speed as a function of the spectral slope a, though, 
strictly speaking, it would be a function of a. 

In our case, we see that, there are two sources of radiation from the disc, (i) the CENBOL 
and (ii) the Keplerian disc. Let us first concentrate on the CENBOL contribution. From Figs. 
(3a) and (3c), we see that for z— > few xlOr^, {£c -\-Vc)^'^J^c, hence Eq. (8) determines 
Veq^l at such distauces. 

In general, in the hard state also, Lk is not negligible. Though at 2;^large, Ek^Fk^Pk, 
but at zr^ few xlOr^, {Sk + Vk) > '^^k- Hence for higher Keplerian luminosity, Veg^l 
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Figure 4. Variation of (a) Veq and (b) ^ with log(z) for I = 0.3 and Xs = 20rj. Various curves correspond to fnx = 0.01 
(solid), rhif = 1 (dashed) and niK = 6 (long-dashed). Variation of (c) Veq and (d) ^ with log{z). for £ = 0.12 and m_fc = 0.5. 
Various curves corresponds Xg = lOr-g (solid), Xs = 20rg (dashed) and Xg = 30rg (long-dashed). 

is achieved at distances around a thousand Schwarzschild radii. In Figs. (4a) and (4b), 
we have plotted Vgq and ^ with log{z). Various curves correspond to rriK ~ 0.01 (sohd), 
niK — l(dashed) and = 6 (long-dashed). We choose I — 0.3 and Xg — 20rg for all 
the plots. We see that close to the black hole, Veq is independent of Ik-, because at such 
distances radiation from the CENBOL dominates. If Lq^Lk {i.e., solid curve), then we see 
that at around lOOr^. For niK = 1 (=^x~0.03, i.e., dashed curve), and for rrix = 6 

{i.e., long-dashed curve), v^q-^l at distances over lOOOr^. One should note that at 2; < l.SSr^, 
Veq < 0. The reason for this can be clearly seen from Fig. (3c), which shows that < 
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at the same distance from the black hole. This means that very close to the black hole 
not only gravity pushes in matter but also radiation would also push matter inside, thus 
vindicating Eq. (10). To clarify this point we have plotted ^ with log{z) in Fig. (4b). Curve 
styles match those of the corresponding case in Fig. (4a) . We see that within the funnel like 
region of the CENBOL, the radiation is completely dominated by the CENBOL itself. And 
also due to this particular geometry ^ > 1, resulting Veq much lesser than the velocity of 
light. In case m^^— ^small (solid), for z > lOOr^, .^^1, resulting Veq^l- In case rriK^ large 
(dashed and long-dashed curves), only at around a thousand Schwarzschild radii, This 
means that, if one increases m^, then the higher velocities achieved due to the acceleration 
of jets by CENBOL photons, might be decelerated. It is quite obvious though, that the 
effect of Keplerian radiation is quite marginal. Of course one should keep in mind that if 
the radiative moments due to CENBOL and Keplerian radiation are comparable at infinite 
distances, then there is a possibility of increased terminal speed, with the increase in mx- 

We have seen from Eq. (9), that ^ not only depends on mK/C- but also on Xs- We now 
investigate the dependence of Veq on Xs- In Fig. (4c), v^q is plotted with log{z) for £ = 0.12 
and triK — 0.5. Various curves represent Xg — lOr^ (solid), Xg — 20rg (dashed) and Xg — SOr^ 
(long-dashed). With the increase in Xg, as the size of CENBOL increases, hence it can only 
behave like a point source (for which ^^1) farther out, and we see that v^q^l farther 
outwards from the black hole. To clarify, we have also plotted ^ with log{z) for i = 0.12 and 
TJiK = 0.5 in Fig. (4d). Various curves represent Xg = lOr^, (solid), Xg = 20rg (dashed) and 
Xg — 30rg (long-dashed). As has been just explained, we see that ^^1 for larger value of z, 
increased. 

3.2 Velocity profile 

So far, we have only discussed about the equilibrium velocity and its dependence on i, tjik 
and Xg. This was done to study the upper limit of allowed velocity as a function of z. We 
have seen that, if CENBOL radiation dominates over that from the Keplerian disc, then jets 
can be accelerated to very high velocities within around lOOr^, but Veq is only a measure of 
velocity which signifies only the domain of radiative deceleration or acceleration, and is not 
the actual velocity. From Eq. (6) it is quite clear that the radiative acceleration and hence 
V itself, depends on all three quantities £, J-', V and not only just on ^. Equation (6) also 
shows, radiative acceleration also depends on v in a very complicated way. To compute actual 
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Figure 5. (a) Variation of v with log{z), for rriK = 0.5. Various curves correspond to f. = 0.36 (solid), i = 0.24 (dashed) and 
£ = 0.12 (long-dashed), (b) Variation of i? with £. Various curves correspond to m/f = 0.4 (solid), rfiK = 1-4 (dashed) and 
tHk = 2.4 (long-dashed). The shock location Xg = 20rg, for both the cases. Injection parameters Zin = 5rg, Vi„ = 10~^. 

velocity, one has to integrate Eq. (6). Jets are believed to be produced from the post shock 
region (Chakrabarti 1999, Das & Chakrabarti 1999, Das et al. 2001). If they are generated 
from the post-shock region then they should start with very small velocity. In this paper we 
are not concentrating on the generation of jets. Thus we just put the injection height to be 
close to the black hole and the injection velocity to be small. We choose injection parameters 
to be Zin = 5rg and f j„ = 10~^. Before discussing the results, let us define terminal speed. 
The terminal speed {-d) is the constant velocity at infinite distances or -d — v\z-*oo- In Fig- 
(5a), three velocity v is plotted with log{z), where ttik — 0.5 and Xg — SOr^,. Various curves 
correspond to ^ = 0.36 (solid), £ — 0.24 (dashed) and i — 0.12 (long-dashed). We see that 
close to the black hole as v<^l, T dominates resulting in a steep rise in v, as v increases 
the jet starts to feel the drag force, resulting in a somewhat less steep increase in v. At 
z > lOOvg, the radiative moments tend to become weak and the jet settles to a constant 
velocity at large z or the terminal speed d. It is clear that, v increases with £ and in the three 
cases depicted, the terminal velocity of jets are ■j?~0.93 (sohd), i?~0.91 (dashed) and ■j?~0.88 
(long-dashed). It may seem curious that even if £ is increased by equal interval the terminal 
speed achieved, does not increase by equal interval. The reason is two fold, first of all from 
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Eq. (6), one can see that the gradient of f is a non-hnear function of the moments. The 
second reason being the existence of 7'^ term in Eq. (6), which suppresses the acceleration, 
as V increases to values close to that of hght. 

In Fig. (5b), the terminal speed -& is plotted with £, for Xg — SOr^,. Various curves 
correspond to rhx — 0.4 (sohd), ttik — 1-4 (dashed) and ttik — 2.4 (long-dashed). The 
terminal speed i? increases with £, though for higher £, the increase of i? decreases, for the 
same two reasons discussed for the previous Figure. We also see that the dependence of 
on tuk is marginal. In Fig. (3), we have seen that the CENBOL contribution to the various 
space dependent part of the radiative moments, is a few orders of magnitude higher than that 
from the Keplerian contribution, hence the marginal dependence of with ttik is expected. 
It is to be noted though, that for £ > 0.15, ■»? decreases with increasing ttik and for £ < 0.15, 
1? increases with increasing tuk, albeit the dependence is very weak. For the higher values 
of i, the jet, powered by the CENBOL radiation, achieves very high velocity within a few 
tens of Schwarzschild radii. But for the next thousand of Schwarzschild radii or so, ^ > 1, 
as {£k + Vk) > '^J^K in that region [see. Fig. (3c-3d)]. When v becomes high, the drag 
force increases, and slows down the jet. At a larger distance, a slight increase in radiative 
moments, due to the increase in triK-i is not sufficient to increase When £ is small, the 
velocity achieved within few tens of Schwarzschild radii is not that high, hence the drag 
force is less. But at large distances, the radiative moments due to CENBOL and Keplerian 
disc becomes comparable hence there is a slight increase in i? with iriK- 

In Fig. (6), 'd is plotted with rriK-i for Xg — 20rg. Various curves correspond io £ — 0.3 
(solid), £ = 0.24 (dashed), £ = 0.18 (long dashed), and £ = 0.12 (dashed-dotted) . We see 
that, as in Fig. (5b), for higher £, ■& decreases with increasing (see, sohd, dashed and 
long dashed curves), but for lower values of £, i? increases with increasing niK (see, dashed- 
dotted curve). So one can conclude that in the hard state, jets arc basically accelerated by 
radiations from CENBOL and the terminal velocity has a weak dependence on radiation 
from Keplerian disc. 

Let us now investigate the dependence on Xg- In Fig. (7a), v is plotted with log{z), for 
£ — 0.18 and ttik — 0.5. Various curves correspond to Xg — lOr^ (solid), Xg — 20rg (dashed) 
and Xg — SOrg (long-dashed). We see that with the increase in Xg, the acceleration is reduced. 
For Xg = lOrg (solid), we see that the jet is experiencing tremendous acceleration and within 
a short distance of z~40rg, it achieves a velocity ~0.9, and gradually settles to a terminal 
value oi "d — 0.93. As the shock location is increased to Xg = 20rg (dashed) and Xg — 30rg 
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Figure 6. Variation of i? with mn, for Xs = 20rg. Various curves correspond to £ = 0.3 (solid), i = 0.24 (dashed) and £ = 0.18 
(long-dashed) and £ = 0.12 (dashed-dotted) . Injection parameters are — ^'^g ■> '^in — 

(long-dashed), acceleration is weaker and the terminal speeds achieved are ■& — 0.9 and 
■d = 0.88 respectively. Thus, the terminal speed depends strongly on shock location and 
decreases with increasing shock location. 

In Fig. (7b), we have plotted d with I. The solid curves correspond to rriK = 1.4 and the 
dashed curves correspond to rriK — 2.4. The shock locations Xg, are marked on each pair 
of sohd and dashed curves {xg — lOr^, 20rp, SOr^) in the Figure. We see that in general ■& 
increases with £. We also see that, for particular values of Xg, rriK has very limited influence 
on ■& (see each pair of sohd and dashed curves marked with values of Xg — 10, 20, 30), and 
also that with increasing x^, itlk is less and less effective in determining d. In particular, 
the solid and the dashed curves marked Xg = 10 are distinguishable, but are increasingly 
less distinguishable for Xg — 20 and Xg — 30. With the increase in Xg, the inner edge of the 
Keplerian disc is increased. We also know that the magnitude of Keplerian disc intensity is 
more, closer to the black hole, i.e., if Xg increases then the size of the Keplerian disc as well 
as Keplerian luminosity decreases. Thus with the increase of Xg, Keplerian radiation will be 
less effective in determining i9. Similar to Fig. (5b), we also notice that if £ > 0.11 for the 
pair of curves (solid and dashed) marked 10 {i.e., = 10), then increases with decreasing 
rriK, but decreases with decreasing ttik, for £ < 0.11. Though this crossing over value of £, 
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Figure 7. (a) Variation of v with log{z) for i = 0.18 and mn = 0.5. Various curves correspond to Xs = lOr^ (solid), Xs = 20rg 
(dashed) and Xs = 30rg (long-dashed), (b) Variation of i? with i. The solid curves correspond to tHk = 1.4 and dashed curves 
correspond to tHk = 2.4. The shock locations Xs, are marked on each pair of solid and dashed curves {xg = Wrg, 20rg, 30rg). 
(c) Variation of i? with Xs, for tHk = 0.5. Various curves represent i = 0.24 (solid), £ = 0.18 (dashed), £ = 0.12 (long dashed) 
and i = 0.06 (dash-dot), (d) Variation of i? with Xs, for i = 0.18. Various curves represent tHk = 1.5 (solid) tHk = 3 (dashed) 
and rhit: = 4.5 (long dashed). Injection parameters Zi„ = 5rg, Vi„ = 10~^. 

increases with increasing Xs, i.e., for Xg = 20rg, the crossing over occurs at £ = 0.15, and for 
Xg = 30rg this occurs at = 0.24. We also see that, generally, decreases with increasing 
Xs- In Fig. (7c), "d is plotted with Xg, for tuk — 0.5. Various curves correspond to £ — 0.24 
(solid), ^ = 0.18 (dashed), £ = 0.12 (long-dashed) and i = 0.06 (dashed-dotted) . It is clear 
that -& decreases with Xg, though for a fixed value oi Xg, ■& increases with I. In Fig. (7d), ■& 
is plotted with x^, for constant values of £ = 0.18. Different curves corresponds to ttik = 1-5 
(solid), rriK = 3 (dashed) and rriK = 4.5 (long-dashed). The first thing to notice is the weak 
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dependence of 'd on ttlk- We also notice that, for Xg < 22rg, i} decreases with increasing rriK 
and for Xg > 22rg, '& increases with increasing tjik- The reason for this is that the CENBOL 
intensity falls with increasing Xs- Thus for a larger Xs the radiative moments due to the 
CENBOL at infinite distances are comparable to that due to the Keplerian disc. Hence with 
increasing rfiK, "& increases, exactly for the same reason as has been presented while discussing 
Fig. (5b). In general, one can conclude from Figs. (7a- 7d), that jets can be accelerated to, up 
and around 90% the velocity of light, provided the x^^lOr^ — SOr^, and £^0.1 — 0.2. As we 
discussed earlier, such values of ^ are not unreasonable when amplification of photon energy 
takes place due to Comptonization and r ~ 1. For smaller r, amplification factor is higher 
(so that £ ~ 1 is achievable) but in that case, both the number of very energy photons goes 
down and also the efficiency to deposite radiative momentum goes down dramatically. 

In the preceding paragraphs of this subsection, we have studied the issue of radiative 
acceleration of jets, and on its dependence on three parameters namely Xs-, d- and rriK- In case 
of CENBOL radiation, the information we have provided is through its total luminosity, but 
Keplerian luminosity is governed by two parameters Xs and rriK [see, Eq. (11)]. To have a 
better understanding, now we study how the relative proportions of CENBOL and Keplerian 
luminosity affect the terminal speed of jet. 

In Fig. (8a), is plotted with Ik/(- — the ratio of Keplerian and CENBOL luminosities, 
for Xa = 20rg. Various curves correspond to i = 0.3 (solid), £ = 0.24 (dashed), i = 0.18 
(long-dashed) and i = 0.12 (dashed-dotted). The ratio of luminosities is kept less than one 
to mimic the hard state of the accretion disc. We see that for i — 0.12, -d increases with the 
increase of Keplerian luminosity. For higher CENBOL luminosities (see, solid, dashed and 
long-dashed curves) we observe that decreases with increasing Keplerian luminosity. This 
has been addressed while discussing Fig. (5b), Fig. (6) and Fig. (7b), i.e., for Xg — 20rg, 
if i? < 0.15, 'd increases with the increase of Keplerian luminosity. It is thus clear that 
if Xs = 20rg then for £ higher than 0.15, jets can be accelerated to very high terminal 
velocities. In order to study the dependence of shock location and the Keplerian luminosity, 
in Fig. (8b), we have plotted -& with i^/i ior i — 0.24. The Keplerian luminosity thus is 
increased from 3% to 75% of CENBOL luminosity. Various curves correspond to Xg — lOvg 
(solid), Xg — 20rg (dashed) and Xg — 3Qrg (long-dashed). For Xg — lOr^, is around 0.93 
but decreases with the increase of Keplerian luminosity. With the increase in shock location 
the terminal velocity is less = 0.91 for Xg = 20rg), but at the same time, decrement of 
the terminal velocity due to the increase of Keplerian luminosity is also less. For Xg — SOrg, 
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'i?~0.895 is lesser, but the change due to the increase of Keplerian luminosity is even less and 
■J? remains almost constant. Thus we conclude that if the shock location is between lOr^ — 
20rg then jets can be accelerated to terminal speeds above 90% of the velocity of light, for 
disc luminosities around 20% of the Eddington luminosity. 



4 DISCUSSION AND CONCLUDING REMARKS 

It is well known that high energy photons can produce particle- antiparticle pairs close to the 
inner edge of a disk. If the photon energy hu ^ 2mc^, then an electron- positron pair may 
be created, where h is the Planck's constant, u is photon frequency and m is the electron 
(or positron) mass. If, on the other hand, electron and positron collide it will annihilate 
each other to produce two Gamma-ray photons, a process called pair annihilation. Clearly, 
to produce electron-positron jets the pair-production process has to dominate over pair 
annihilation. Workers (e.g. Mishra and Melia, 1993; Yamasaki, Takahara and Kusunose, 
1999) discussed the production of electron-positron pairs from inner part of the accretion 
disc. Observationally, there arc reports of pair dominated jets (Sunyaev et al. 1992; Mirabel 
and Rodriguez 1998; Wardle et al. 1998) from galactic black hole candidates to quasars. 
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Though there is httle doubt on the existence of pair dominated jets, radiative acceleration 
of such jets on the other hand is a different issue altogether. If the number density of pairs 
created around the the black hole is too high, then radiative acceleration would be ineffective. 
In this present paper, we have ignored the details of formation of the pair plasma jets and 
have only concentrated on the radiative acceleration of optically thin pair dominated jets. 

In this paper, we supplied the CENBOL intensity (Ic), and the shock location Xs as 
free, but reasonable parameters. We have separately treated two velocity variables, (i) the 
equilibrium velocity Veq and (ii) the actual velocity v. While Veq decides how much velocity 
is allowed before deceleration sets in, v gives us what 'net' value of velocity is achieved 
by actual acceleration. We have shown that if only CENBOL radiation dominates over the 
radiation from the Keplerian disc, i'eg~l is achieved within about few tens of Schwarzschild 
radii, but the same condition is achieved at over a thousand Schwarzschild radii for much 
higher values of Keplerian luminosity. From Figs. (4a-4b), we have seen as the Keplerian 
accretion rate is increased, Vgq decreases in the range z—^15rg — ISOOrg. 

From Fig. (3d), we have seen that {Ek + Pk) > '^Fk in this very range, so in this range 
the radiation drag is higher due to Keplerian radiation. This means that for higher values 
of I and for fixed Xg i.e., for higher CENBOL intensity, the flow will tend to achieve very 
high velocities within few tens of Schwarzschild radii, but at the same time if Keplerian 
luminosity or ^K is increased, then within the range ISr^ — ISOOr^, Keplerian radiation will 
tend to reduce the high velocities because of increased radiation drag. 

From this study we generally conclude that, 

(i) radiative acceleration of electron-positron jets does achieve relativistic terminal speed. 

(ii) The space dependant part of the radiative moments from the post-shock region, domi- 
nates the corresponding moments from the Keplerian disc. 

(iii) In general, the terminal speed of jets increases with increasing post-shock luminosity. 

(iv) Post-shock radiative intensity decreases with increasing shock location, so the terminal 
speed also decreases with increasing shock location. 

(v) Keplerian radiation has marginal effect in determining the terminal speed. 

(vi) Our calculations show that, if the shock in accretion is located at around 10 — 20 
Schwarzschild radii, and if the post-shock luminosity is about 10% to 20% of the Eddington 
luminosity, then electron-positron jets can be accelerated to terminal speeds above 90% the 
speed of light. 
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